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describing  the  program  structure  and  flow  with  complete  Samples 
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PREFACE 


In  April  1973,  the  FACT  model  was  designated  as  the 
Navy  Interim  Standard  Model  for  the  prediction  of  trans¬ 
mission  loss  in  an  environment  which  can  be  characterized 
by  a  flat  bottom  and  a  single  sound  speed  profile.  At 
that  time,  preliminary  documentation  was  distributed  to 
Navy  laboratories  and  other  qualified  users. 

These  volumes  contain  a  detailed  description  of  the 
physics,  mathematics,  and  computer  implementation  of  the 
FACT  program.  Questions  concerning  the  model  or  suggestions 
for  its  improvements  should  be  addressed  to  the  Acoustic 
Environmental  Support  Detachment  of  the  office  of  Naval 
Research. 
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J  2  .  , 
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1.0  INTRODUCTION  AND  BACKGROUND 

On  April  1,  1973,  the  Acoustic  Environmental  Modeling 
Coordinator,  Dr.  J.  B.  Hersey,  designated  a  new  Navy  Interim 
Standard  Transmission  Model  (FACT  -  Fast  Asymptotic  Coherent 
Transmission)  for  ocean  environments  which  may  be  treated 
within  the  context  of  a  single  sound-speed  profile  and  a  flat 
bottom.  This  model  resulted  from  a  three-year  effort  in 
acoustic  modeling  sponsored  by  the  Long  Range  Acoustic  Pro¬ 
pagation  Project  (LRAPP)  within  several  Navy  and  industrial 
laboratories.  As  a  part  of  thin  program  LRAPP  hosted  in 
March  1971  a  workshop  on  acoustic  modeling  by  ray-tracing 
techniques.  Sixteen  different  organizations  active  in  the 
field  of  ray  tracing  participated,  and  a  by-product  of  the 
workshop  was  a  "consensus  model"  consisting  of  components 
drawn  from  several  models.  Subsequent  to  the  workshop  Bell 
Telephone  Laboratori  •  ;  (BTL)  developed  for  LRAPP  a  special 
purpose  program  base*-,  upon  this  consensus  model  and  oriented 
specifically  to  the  Airborne  Sensor  Range  Prediction  (ASRAP) 
program  at  Fleet  Numerical  Weather  Central  (FNWC  ,  Monterey, 
Cal.)  where  speed  was  crucial.  In  October  1972,  a  preliminary 
version  of  this  model  was  turned  over  to  AESD  which  in  turn 
delivered  to  FNWC  on  March  15  a  fully  automated  upgraded  model 
For  a  typical  deep-water  ASRAP  case  the  transmission  loss  is 
computed  every  0.5  nm  to  125  nm  at  four  frequencies  for  each 
pair  of  source  and  receiver  depths  in  4  to  8  seconds. 

While  the  FACT  model  was  originally  designed  for  ASRAP, 
it  is  sufficiently  general  for  many  other  applications. 

Since  turning  it  over  to  FNWC,  AESD  has  distributed  it  to 
37  other  Navy  and  industrial  activities  in  six  different 
countries.  The  documentation  of  this  model  has  to  date 
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cons  in  tod  of  four  HT1.  pub  l  ications  (Spotford,  1972a,  1972b, 
1972c  and  1971a)  and  descript  ive  material  generated  by  AKSD 
and  included  as  part  of  a  fully  annotated  listing.  This 
report  represents  the  final  technical  documentation  of  the 
FACT  model  ns  of  1  December  1973.  As  with  any  computer  re¬ 
presentation  of  a  complex  physics  problem  substantive  mod¬ 
ifications  and  improvements  in  the  physics  are  likely,  as 
are  some  bugs  in  the  code  which  have  eluded  initial  testing. 
The  FACT  model  has  boon  extensively  tested  and  is  run  several 
hundred  times  daily  by  FNWC.  As  bugs  are  uncovered,  however, 
AF.SD  will  provide  correction  sets  to  all  users.  In  addition, 
when  any  improvements  in  physics  are  made,  technical  addenda 
to  this  report  will  be  issued. 

The  FACT  model  is  a  ray-acoustics  model  designed  for  the 
computation  of  transmission  loss  as  a  function  of  range  and 
frequency  .it  fixed  source  and  receiver  depths.  The  classical 
ray  treatment  has  been  augmented  with  higher  order  asymptotic 
corrections  and  the  phased  addition  of  selected  paths.  These 
improvements  yield  a  more  complete  model  of  diffraction  and 
coherence  effects,  respect  i  ve  1  y ,  found  to  be  important,  in 
long-range  Inw-f requeney  propagation.  The  basic  physics 
package  has  been  fully  automated  so  that  the  only  inputs  re¬ 
quired  are  the  specification  of  t:lu>  environment  and  the  para¬ 
meters  of  tlu»  case  desired.  While  programs  incorporating  more 
precise  ray  tracing  may  exist  (involving  curvilinear  sound 
speed  profiles,  ray  targeting,  etc.),  none  appears  to  contain 
as  complete  and  automated  a  treatment  of  the  diffraction  and 
coherence  effects  so  important  at  low  frequencies.  The  pro¬ 
gram  has  been  constructed  in  an  extremely  modular  form  which 
should  accomodate  both  the  incorporation  of  new  physics 
packagos  as  well  as  the  need’s  of  the  R&D  community  for  a 
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flexible  basic  framework  within  which  new  theories  and  tech¬ 
niques  may  be  investigated. 

This  report  consists  of  four  basic  sections  in  two 
volumes : 

Volume  T 


1.  A  summary  of  the  capabilities  of  the  FACT  model, 
its  strengths  and  weaknesses,  and  those  areas  in  which 
future  improvements  are  planned  as  part  of  a  detailed  des¬ 
cription  of  the  physics  represented  by  the  program.  Ex¬ 
amples  are  provided  comparing  normal-mode  computations  and 
FACT  output  to  i.1  'us urate  several  of  its  novel  features. 

2.  A  description  of  the  mathematical  implementation  of 
the  model  with  all  relevant  egaations. 

Volume  IX  (Provided  routinely  to  all  recipients  of  the 
FACT  Program,  available  upon  request,  and  taken  nearly 
verbatim  from  a  report  issued  by  Ocean  Data  Systems,  Inc., 
documenting  the  FACT  program  (Da^er,  1973) 1 

1.  A  description  of  the  program  structure  and  the  pro¬ 
gram  flow  including  readable  flow  cnarts  of  all  routines. 

2.  A  copy  of  the  original  FACT  documentation  prepared 
by  AESD  and  normally  included  with  a  full  program  listing 
("FACT  Handout").  Thi3  documentation  included  samples  of 
the  complete  input  and  output  for  three  teat  cases  to  assist 
new  users  in  implementing  and  debugging  the  program  on  thair 
computers.  The  full  program  listing  hat  not  been  included 


in  this  report.  Upon  n-g...-yi  .  -  Ai...i  .  ...  i  PACT  listing 
and  a  punched-card  deck  arc  av.ii’' 
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Finally,  the  author  is  especially  indebted  to  Mr.  c. 
L.  Baker  of  Ocean  Data  Systems,  Incorporated,  who  upgraded 
the  original  computer  code  to  a  product  with  considerably 
fewer  bugs  and  in  a  form  which  is  sufficiently  modular  to 
be  useful  to  a  wide  community.  Mr.  Baker  converted  a  pro¬ 
gram  written  by  scientists  and  engineers  into  a  computer 
program,  and  is  entirely  responsible  for  the  more  lucid 
portions  of  the  FACT  code. 


2.0  SUMMARY  OF  BASIC  CAPABILITIES  OF  THE  FACT  MODEL 


This  section  is  intended  to  familiarize  the  reader  with 
the  basic  capabilities  and  limitations  of  FACT,  in  terms  of 
the  physical  model  it  represents.  A  discussion  is  also  in¬ 
cluded  of  those  areas  of  the  program  which  will  be  improved 
and  extended  in  the  future.  It  should  be  re-emphasized 
that  the  modularity  of  the  program  will  permit  modelers  to 
test  new  theories  and  approaches  quite  easily  using  the  basic 
FACT  framework.  The  discussion  in  this  section  has  been 
divided  into  two  major  parts  addressing: 

1.  The  model  used  for  the  ocean  environment; 


2.  The  ray  treatment  of  this  environment. 
2.1  The  Ocean  Environment 


The  model  used  for  the  ocean  environment  is  "layered"; 
that  is,  neither  the  sound-speed  profile  nor  the  bottom 
depth  are  range  dependent.  The  sound-speed  profile,  c(z)  , 
is  approximated  by  a  continuous  piecewise-linear  function 
of  depth  z,  so  that  within  each  layer  the  sound-speed  gra¬ 
dient,  c*  is  a  constant  and  at  layer  boundaries  c'  is 

discontinuous.  (The  effects  of  such  discontinuities  have 

been  well  studied  (Pedersen,  1961)  and  will  be  addressed  in 
sections  3.1  and  3.2.)  Within  the  water  column  a  frequency- 
dependent  volume  attenuation  (or  absorption)  is  included 
of  the  form 
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where  f  is  the  frequency  in  kilohertz  and  o  is  a  loss  in 
dB  per  nautical  mile.  This  expression  corresponds  to 
Thorp's  absorption  (Thorp,  1965}  above  1  kHz  and  is  an 
empirical  result  derived  by  FNWC  (Parka  I,  1969)  below 
1  kHz. 

The  ocean  bottom  is  modeled  as  a  reflecting  surface 
where  each  ray  suffers  a  reflection  loss  which  is  dependent 
upon  grazing  .ingle,  frequency,  and  bottom  class.  The  loss 
suffered  as  a  function  of  these  parameters  is  given  by  the 
FNWC  bottom-loss  tables  (Bassett  and  Wolff,  1970)  for 
frequencies  loss  than  1  kHz  or  greater  than  3.1  kHz,  and  by 
the  new  Naval  Oceanographic  Office  (NOO)  Navy  Standard 
curves  from  1  kHz  to  3.5  kHz  (Christensen,  Frank  and  Kaufman, 
1972)  . 


The  ocean  surface  is  modeled  as  a  perfect  reflector 
'with  a  180-degree  phase  shift)  for  all  rays  considered. 

For  propagation  in  a  surface  duct  rays  are  not  used  and  the 
surface-duct  module  (to  be  described  subsequently )  include? 
a  rough-surface  loss  which  is  dependent  on  frequency,  sea 
state  and  mixed  layer  depth  (i.e.  duct  thickness). 

In  the  three  areas  of  volume  absorption,  bottom  reflec¬ 
tivity,  and  rough-surface  loss  improvements  are  expected. 

The  volume  absorption  is  currently  added  to  the  total  trans¬ 
mission  loss  as  a  function  of  horizontal  range.  Ideally 
volume  absorption  should  be  computed  on  a  ray-by-ray  basis 
using  the  arc  length  of  the  ray  for  the  appropriate  distance. 
It  is  expected  that  this  more  accurate  computation  may  be 
available  in  the  near  future,  though  it  will  increase  the 
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running  time  somewhat.  More  significantly,  the  low- 
frequency  absorption  will  be  updated  as  soon  as  the  R&D 
community  agrees  on  appropriate  values  applicable  to 
individual  ray  paths.  A  more  consistent  merger  of  the  FNWC 
and  NOO  bottom- Loss  types  is  needed  as  well  as  considerably 
more  data  on  low-frequency  reflectivity.  Similarly  a  good 
surface-reflectivity  model  (dependent  presumably  on  grazing 
angle,  frequency,  and  some  measure  of  surface  roughness)  is 
needed,  again  directly  applicable  to  the  coherent,  specular 
scatter  of  individual  rays. 

2.2  The  Basic  FACT  Model 


This  section  contains  a  qualitative  description  of  the 
basic  physics  in  the  FACT  model  -  that  is:  what  effects  it 
does  and  does  not  treat;  what  the  relative  levels  of  accuracy 
of  the  various  approximations  employed  are;  and  where  treat¬ 
ments  might  be  extended  or  refined  in  the  future.  The  de¬ 
tailed  discussion  of  how  these  computations  are  performed 
(equations,  approximations,  algorithms,  etc.)  are  reserved 
for  Section  3. 

Acoustic  propagation  in  the  ocean  is  governed  by  the 
wave  equation.  Ray  tracing  is  a  convenient  and  powerful 
technique  for  approximating  the  exact  wave  solution.  The 
traditional  ray  solutions  determine  all  relevant  ray  paths 
connecting  the  source  and  receiver  and  sum  the  geometric 
intensities  of  these  paths  to  obtain  the  total  transmission 
loss.  This  traditional  approach  is  already  an  approximation 
to  riqorous  ray  theory  which  requires  the  phased  sum  of  the 
ray  amplitudes.  The  justification  for  this  approximation  is 
that  the  phased  sum  oscillates  quite  rapidly  with  range,  and 
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that  uncertainties  in  the  geometry  and  environment  preclude 
an  accurate  estimate  of  phase.  Hence  the  phased  sum  is  re¬ 
placed  by  the  rms  sum  as  an  average  over  the  uncertainities 
in  position  and  environment.  As  will  be  shown  subsequently, 
in  many  significant  cases  this  rms  sum  is  not  representative 
of  these  or  any  other  meaningful  averages. 

A  second  area  where  traditional  ray  tracing  requires 
improvements  is  in  the  computation  of  individual  ray  inten¬ 
sities.  According  to  classical  ray  theory  the  intensity  of 
a  ray  (or,  more  appropriately,  a  ray  bundle)  is  inversely 
proportional  to  the  geometric  spreading  or  divergence  of  the 
bundle.  In  the  ocean  environment  focal  points  and  surfaces 
occur  where  the  ray  divergence  is  zero  and  infinite  ray 
intensities  are  predicted,  usually  accompanied  by  immediately 
adjacent  areas  of  zero  intensity  (i.e.  no  rays).  The  acous¬ 
tic  field  associated  with  these  focusing  regions,  or  caustics, 
is  improperly  modeled  by  classical  ray  theory.  The  error  at 
low  frequencies  may  be  substantial  over  large  regions. 

Recent  advances  in  ray  theory  and  the  analysis  of 
acoustic  measurements  have  permitted  these  two  major  defic¬ 
iencies  of  classical  ray  theory  to  be  overcome.  The  FACT- 
model  treatment  of  these  improvements  as  well  as  algorithms 
for  approximating  low— frequency  cut-off  effects  and  axis- 
to-axis  propagation  are  described  in  the  following  subsections. 

2.2.1  Coherent  Summation  of  Paths 

At  even  very  low  frequencies,  the  coherent  summation  of 
all  ray  paths  produces  very  rapid  oscillations  in  the 


transmission  loss  as  a  function  of  range.  For  most  applica¬ 
tions  these  oscillations  are  too  rapid  to  be  useful  in  detail 
and  should  be  smoothed.  Also  uncertainties  in  the  environ¬ 
ment  are  sufficient  to  preclude  the  prediction  of  the  precise 
position  of  these  peaks  and  nulls.  What  is  needed  than  is 
a  range-averaging  technique  which  smooths  these  rapid  fluct¬ 
uations  while  preserving  long-term  significant  departures 
from  the  rms  sum  of  all  paths. 

The  long-term  departures  from  the  rms  sum  are  primarily 
due  to  long-range  surface-image  interference  effects.  Fig¬ 
ure  2-1  illustrates  a  typical  pair  of  RSR  paths  which  differ 
only  in  the  vicinity  of  the  source.  The  phase  difference 
between  these  paths  consists  of  the  travel-time  difference 
and  a  phase  reversal  of  n  for  the  path  with  the  additional 
sui  t  ace  reflection.  For  shallow  sources  and  low  frequencies 
this  phase  difference  changes  quite  slowly  with  range  and 
may,  over  entire  convergence  zones,  change  by  considerably 
less  than  1  full  cycle  (DeAngelis  and  Spofford,  1970  and 
Hanna,  1970) ,  Hence  the  phased  sum  of  these  two  paths  may 
potentially  differ  by  a  uniform  amount  from  the  rms  sum  over 
significant  range  intervals. 


Fortunately  the  phase  difference  may  be  estimated  quite 
accurately  from  the  ray  geometry  in  the  immediate  vicinity 
of  the  source  and  detailed  travel-time  calculations  are  not 
necessary.  Since  the  two  ray  amplitudes  are  essentially 
equal  a  rapidly  computed,  local  phaseu  sum  may  be  obtained. 
This  approach  has  been  incorporated  into  the  FACT  model. 
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Figures  2-2,  2-3,  2-4  and  2-5  compare  predictions  for 
an  RSR  convergence  zone  using  FACT  ("RAY  THEORY")  with  an 
exact  normal-mode  result  for  frequencies  of  25,  50,  100  and 
200  Hz  respectively.  (The  leading  edge  of  this  convergence 
zone  contains  a  caustic  which  will  be  discussed  subsequently.) 
The  shallow  source  experiences  very  strong  surface-image 
interference.  As  the  frequency  increases  the  net  effect 
shifts  from  uniform  near  cancellation  at  25  Hz  to  uniform 
enhancement  at  50  and  100  Hz.  At  200  Hz  the  first  null  in 
the  Llbyd's  Mirror  interference  pattern  occurs  at  33  miles 
and  is  predicted  fairly  accurately  by  FACT.  It  should  be 
noted  that  the  more  rapid  oscillations  in  these  mode  pre¬ 
dictions  are  being  intentionally  smoothed  by  FACT. 

In  addition  to  either  this  fully  coherent  combination 
of  the  two  paths  or  the  incoherent  rms  sum,  a  third  "semi- 
coherent"  option  is  available.  The  FACT  program  computes 
transmission  loss  at  discrete  range  points,  and  for  deeper 
sources  and/or  higher  frequencies  than  those  of  Figures 
2-2,  2-3,  2-4  and  2-5,  the  interference  pattern  may  be 
undersampled  yielding  an  apparently  random  structure.  Since 
the  range  step  is  indicative  of  the  amount  of  detail  usually 
desired  in  a  prediction,  the  semi-coherent  option  provides 
an  automated,  smooth  transition  from  the  fully  coherent  two- 
path  svim  to  the  incoherent  sum  as  the  number  of  range  points 
per  cycle  of  interference  decreases  from  6  to  2J/j.  This 
option  is  recommended  for  general  use,  keeping  in  mind,  how¬ 
ever,  that  it  may  yield  range-step  dependent  transmission 
loss.  Figure  2-6  illustrates  the  smoothing  process  in  the 
first  bottom-bounce  region  for  varying  sampling  rates. 

Finally,  these  coherence  effects  are  possible  at  both  source 
and  receiver  if  the  receiver  is  sufficiently  shallow,  and 
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where  appropriate  are  included  at  one  or  both  ends  of  the 
ray  path. 

2.2.2  Caustic  Processing 

There  are  three  types  of  focusing  regions  which  are 
treated  by  FACT:  smooth  caustics,  cusped  caustics,  and 
combined  srnooth/cusped  caustics.  The  treatments  of  the  first 
two  are  based  upon  thorough  analyses  of  the  ray  geometries, 
and  the  errors  associated  with  any  approximations  are  known. 
The  third  case  has  not  been  thoroughly  a.nalyzed,  and  the 
errors  associated  with  its  treatment  are  less  well  defined. 

2.2.2. 1  Smooth  Caustics  (Two- Ray  Systems) 

Figure  2-7  depicts  a  family  of  rays  in  a  simple  pressure- 
gradient  profile  which  form  a  smooth  caustic  after  the  first 
surface  reflection.  The  caustic  is  the  envelope  of  this 
family  of  rays  and  forms  the  classical  boundary  between  the 
illuminated  region  (with  two  rays  through  each  point)  and  the 
shadow  zone  (witn  no  rays  present) .  The  classical  field 
increases  in  the  illuminated  region  to  infinity  on  the  caus¬ 
tic  surface  and  is  zero  m  the  shadow  zone.  The  actual 
field  as  given  by  wave  theory  reaches  a  finite  maximum  value 
just  inside  the  illuminated  region  and  decays  ex|x>nent  ial  ly 
into  the  shadow  zone.  The  key  questions  lets,  unresolved  by 
classical  ray  theory  are  the  values  of  this  peak  amplitude 
and  the  rate  of  decay  (which  at  low  frequencies  may  be  quite 
g  radua l ) . 

Higher  order  asymptotic  estimates  for  these  quantities 
may  be  obtained  (Brekhovskikh,  1960a,  Sachs  and  Si  Unger. 


1971  arul  Ludwig,  1966)  from  the  local  ray  geometry  near  the 
caustic.  FACT  computes  these  expressions  and  yields  the 
shadow- zone  field,  the  peak  field,  and  a  range  average  of 
the  very  rapid  two-ray  interference  in  the  illuminated  region. 
The  convergence  zone  comparisons  of  Figures  2-2,  2-3,  2-4 
and  2-5  correspond  to  the  geometry  of  Figure  2-8,  where  the 
receiver  ("SOURCE"  invoking  the  usual  reciprocity  assump¬ 
tions)  is  deep  and  the  source  is  located  50  feet  below  the 
surface.  At  the  source  depth  the  caustics  of  the  direct 
and  surface-reflected  rays  combine  on  a  phased  basis  and 
the  resulting  sur face- image  interference  effects  must  be 
included.  The  primary  effect  of  the  caustic  corrections  in 
Figures  2-2,  2-3,  2-4  and  2-5  is  seen  in  the  increase  in 
slope  of  the  leading  edge  of  the  zone  as  the  frequency  in- 
eroases.  The  changes  in  peak  level  with  frequency  are  domin¬ 
ated  by  coherence  effects  in  this  case  and  do  not  reflect 
the  weaker  frequency  dependence  of  the  caustic  itself. 

At  very  low  frequencies,  or  for  caustics  which  cross 
the  receiver  depth  nearly  horizontally,*  the  shadow-zone 
fall-off  is  quite  slow  and,  in  the  approximation  of  the  local 
caustic  geometry,  may  be  computed  to  extend  to  great  range. 
Smooth  caustics  typically  originate  or  terminiate  in  a  cusp 
(point  C  in  Figure  2-8  which  will  be  discussed  in  more  de¬ 
tail  in  the  following  section) .  Since  the  shadow-zone 
field  essentially  extends  no  farther  than  the  cusp,  FACT 
tapers  the  field  to  be  fully  decayed  at  the  range  of  the  cusp, 
preserving  its  initial  behavior  at  the  smooth  caustic. 


*  FACT  is  unable  to  include  the  effects  of  fields  assoc¬ 
iated  with  caustics  which  do  not  reach  the  receiver  (or 
ource)  -  for  example  leakage  into  or  out  of  a  surface  duct. 
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2. 2.2.2  Cusped  Caustics  (Three- Ray  Systems) 

Just  as  the  geometrical  acoustics  amplitude  is  invalid 
at  smooth  caustics,  the  smooth-caustic  amplitude  breaks 
down  at  cusped  caustics.  The  cusped  caustic  C,  illustrated 
in  Figure  2-8  occurs  where  two  smooth  caustics  terminate. 
Inside  the  cusp  there  are  three  rays  through  any  point  and, 
as  with  smooth  caustics,  the  field  near  the  cusp  may  be 
determined  from  the  local  ray  geometry  (Holford,  1972). 

An  automated  computation  of  the  field  for  general  cusped 
caustics  is  quite  involved;  however,  the  observation  that 
in  layered  media  a  cusp  always  forms  at  the  same  depth  as 
the  source  simplifies  the  calculation  significantly.  This 
case  is  of  particular  interest  for  shallow-source,  shallow- 
receiver  geometries  (e.g.  sonobuoys)  and  has  been  accounted 
for  in  FACT  (averaging  over  the  rapid  three-ray  interference 
inside  the  cusp  as  in  the  smooth-caustic  interference  region) 

Figures  2-9,  2-10  and  2-11  compare  FACT  ("RAY  TRACING" ) 
with  a  normal-mode  result  for  the  geometry  of  Figure  2-8  at 
25,  50  and  100  Hz  respectively.  The  first  cusp  occurs  at 
~21  nm  and  repeats  every  21  nm.  Each  cusp  is  preceded  at 
35,  58  and  80  nm  by  a  smooth  caustic  which  originates  at  the 
preceding  cusp  and  is  reflected  from  the  surface.  For  these 
convergence  zones  the  smooth  and  cusped  caustics  are  part  of 
the  same  ray  aperture  and  a  potential  four-ray  system  exists. 
The  treatment  of  such  cases  is  described  in  the  next  section. 


2. 2. 2. 3  Smooth  and  Cusped  Caustics  (Four-Ray  Systems) 


Figure  2-12  illustrates  the  canonical  cusped-caus tic 
geometry  foi  a  shallow  source  in  a  warm-water  profile.  For 


a  receivin'  .it  the  source  depth  the  cusp,  C,  is  preceded 
immediately  by  a  smooth  caustic  (emanating  from  the  cusp 
at  C').  for  more  typical,  stronger  thermocline  gradients 
the  separation  is  considerably  less.  Mso  for  deeper 
water  than  shown  in  the  figure  the  up-going  rays  asso¬ 
ciated  with  the  smooth  caustic  would  penetrate  the  cusp 
resulting  in  four  rays  through  points  within  the  cusp. 

In  principle,  this  is  a  basic  four-ray  system  for 
which  the  necessary  asymptotic  expressions  are  not  available. 
In  practice,  the  smooth-  and  cuspcd-caustic  fields  may  be 
superimposed  to  yield  reasonable  results.  For  well- 
separated  systems  as  in  Figures  2-9,  2-10  and  2-11  the 
smooth-  and  cusped-caustic  fields  may  be  added  on  an  rms 
basis.  This  is  the  technique  currently  incorporated  in  FACT. 
For  very  tight  geometries  a  phased  sum  of  the  two  fields  has 
been  required  on  occasion  (Holford  and  Spofford,  1972); 
however,  this  is  a  difficult  computation  to  automate. 

Figures  2- l 3,  2-14  and  2-15  compare  FACT  ("RAY  TRACING") 
and  normal-mode  predictions  for  a  simplified  Pacific  profile 
at  frequencies  of  25,  50  and  100  Hz  respectively.  The  ray 
geometry  is  similar  to  that  in  Figure  2-12,  however,  the 
cusped  and  smooth  caustics  are  much  closer.  Only  the  paths 
corresponding  to  the  RR  rays  have  been  included  in  the  com¬ 
putations  except  for  the  FACT  predictions  from  5  to  30  miles 
which  include  bottom- ref iected  paths.  The  frequency  depen¬ 
dence  of  the  agreement  between  these  predictions  suggests 
that  a  fully  phased  combination  might  yield  better  agreement. 
This  modification  is  currently  under  investigation. 

2.2.3  Low-Frequency  Cut-Off  Effects 

In  the  FACT  model  the  total  transmission  loss  is  com¬ 
puted  by  summing  on  an  incoherent  or  rms  basis  the  intensities 
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of  sets  of  rays  which  contribute  to  the  field  at  each 
range/depth  point  of  interest.  Within  each  set  (to  be  des¬ 
cribed  subsequently)  certain  coherent  combinations  of  paths 
may  have  already  been  performed  either  explicitly  (for  sur¬ 
face-image  interference)  or  implicitly  (near  caustics) . 

The  subsequent  incoherent  combination  of  these  sets  assumes 
that  the  relative  phase  differences  between  sets  are  both 
unpredictable  and  rapidly  changing  with  range.  For  very  low 
frequencies  (geometries  with  dimensions  of  several  wave¬ 
lengths)  both  of  these  assumptions  may  be  incorrect.  Most 
importantly,  as  frequencies  decrease  to  near  cut-off  for  the 
first  trapped  mode,  large-scale  cancellations  occur  resulting 
in  significant  uniform  degradations  from  the  rms  intensity. 

In  the  FACT  model  this  effect  is  approximated  by  re¬ 
ducing  the  amplitudes  of  rays  which  would  experience  uni¬ 
form  destructive  interference.  The  rays  expected  to  exper¬ 
ience  this  interference  are  those  with  angles  shallower  than 
the  ray-equivalent  of  the  first  propagating  mode  as  deter¬ 
mined  by  the  standard  WKB  approximations  (Brekhovskikh ,  1960b). 
This  approach  is  admittedly  approximate  and  attempts  only  to 
capture  very  gross  features.  The  extension  of  ray  theory 
to  a  situation  so  clearly  in  the  domain  of  wave  techniques 
is  speculative  at  best.  Further  study  is  underway;  however, 
until  a  thoroughly  substantiated  technique  is  developed, 
wave  programs  should  be  used  wherever  possible  for  these 
cases. 

2.2.4  Axis-to-Axis  Propagation 

A  persistent  problem  in  ray  tracing  has  been  modeling 
so-called  "axis- to-axis"  transmission.  The  difficulty  is 
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illustrated  in  Figure  2—16  whore  rays  are  traced  from  a 
source  located  on  the  axis  of  a  bilinear  profile.  At  any 
range  on  the  axis  there  exists  an  infinite  number  (as  the 
ray  angle  approaches  zero)  of  refracted  rays  (each  of  non¬ 
zero  intensity)  connecting  the  source  and  receiver.  Hence, 
the  rms  summation  of  these  paths  yields  an’ infinite  inten¬ 
sity  and  the  ray  solution  is  invalid.  Ignoring  low-frequency 
cut-off  effects,  this  problem  is  an  undesirable  by-product 
of  the  linearly  segmented  velocity  profile.  For  profiles 
with  a  continuous  first  derivative  of  e(z)  at  the  axis  (and 
non-zero  second  derivative)-  the  near-axial  rays  focus  period¬ 
ically  at  cuspcd  caustics  along  the  axis  where  the  focal  dis¬ 
tance  is  determined  by  the  second  derivative  of  c(z). 

The  linearly  segmented  profile  is  an  approximation  to 
the  smoother  profile  which  is  usually  modeled  as  having  a 
discontinuous  second  derivative  across  the  axis.  This  model 
simulates  the  basic  asymmetry  of  the  sound-channel  axis 
where  the  decrease  in  sound  speed  associated  with  the 
the rmoc line  gradient  is  overtaken  by  the  increase  in  sound 
speed  associated  with  the  "pressure  gradient".  The  discon¬ 
tinuous  second  derivative  gives  two  weak  focal  points  along 
the  axis  (F,  and  F2  in  Figure  2-17)  between  each  strong  focal 
point  (FQ).  All  rays  about  the  axis  are  focussed  at  FQ, 
whereas  only  upgoing  and  only  downgoing  rays  are  focused  at 
F,  and  F2  respectively. 

The  motivation  for  this  model  for  the  sound-channel  axis 
is  to  capture  the  character  of  transmission  associated  with 
the  near-axial  rays.  The  prominent  feature  of  these  rays  is 
the  strong  focus  FQ.  The  implementation  of  this  approach 
in  FACT  is  accomplished  by : t 


1.  Estimating  the  second  derivatives  of  c(z)  above 
and  below  the  axis  and  establishing  an  equivalent  smooth 
profile; 

2.  Computing  the  corresponding  period  (distance  to 
F0)  for  the  axial  ray  in  the  smooth  profile; 

3.  Finding  the  ray  of  minimum  angle  in  the  1  inear ly 
segmented  profile  with  this  period; 

4.  If  source  and  receiver  are  between  the  upper  and 
lower  turning  points  of  this  ray,  moving  both  the  source 
and  receiver  to  the  depth  of  either  the  upper  or  lower 
turning  point. 

The  resulting  transmission  loss  for  this  set  of  near- 
axial  rays  will  then  exhibit  the  strong  focusing  of  a  cus- 
ped  caustic  (simulating  the  focusing  at  F0)  at  the  appro¬ 
priate  range.  By  moving  the  source  and  receiver  to  the 
same  depth  whenever  they  are  both  between  the  turning 
points,  the  gross  features  of  the  near-axial  wave  field  are 
assumed  to  be  essentially  independent  of  depth  for  a  near¬ 
axis  source. 

2.2.5  Mode l l nq  Surface-Ducted  Propagat ion 

FACT  presently  contains  the  same  surface-duct  module 
as  the  Navy  Standard  Bong-Ranee  Transmission  Model,  RP?0. 
The  model  was  developed  f^i  F\ivC  by  Clay  and  Tatro  (Clay, 
1968)  to  characterize  the  gross  features  of  ducted  propaga¬ 
tion.  The  intensity  in  the  surface  duct  is  found  from  con¬ 
servation  oF  energy  modified  by  additional  losses  (propor¬ 
tional  to  range)  caused  by  duct  leakage  and  rough-surface 
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scattering  of  energy  from  the  duct  (Marsh  and  Schulkin, 

1967).  This  basic  intensity  is  independent  of  source  or 
receiver  depth  so  long  as  both  are  in  the  duct.  For 
cross-duct  cases  (only  one  in  the  duct)  the  intensity  is 
reduced  by  10  dB,  and  for  neither  in  the  duct  no  ducted 
contribution  is  computed. 

A  preliminary  analysis  of  this  model  (Labianca,  1972) 
has  shown  that  the  leakage  term  does  not  correspond  to  the 
mechanism  which  dominates  duct  leakage  and,  in  particular, 
has  the  wrong  dependence  on  the  below-layer  ( thermocline) 
gradient.  An  intensive  investigation  of  all  available  sur¬ 
face  duct  models  has  been  completed  recently  by  NUC  for 
LRAPP  (Morris,  1974) .  This  study  compared  several  models 
with  each  other  and  the  SUDS  transmission-loss  data  (Cummins, 
1972).  Since  none  of  the  candidate  models  emerged  as 
clearly  superior,  AESD  is  currently  studying  possible  mod¬ 
ifications  to  the  surface-duct  equations. 

In  the  normal  operating  mode  no  rays  are  traced  in  the 
surface  duct  and  the  intensity  associated  with  ducted  paths 
xs  computed  by  the  duct  module.  The  user  has  the  option  of 
tracing  rays  in  the  duct  and  computing  intensities  using 
the  standard  FACT  algorithms.  This  procedure  is  not  recom¬ 
mended  in  general  for  several  reasons: 

1.  Long-range  computations  require  excessive  computer 

time; 

2.  No  leakage,  or  surface-scattering  effects  are  in¬ 
cluded; 


3.  Cross-layer  coupling  is  not  computed. 

Nevertheless  comparisons  of  rays  versus  mode  calcula¬ 
tions  can  be  informative  in  cases  where  the  above  considers 
tions  are  not  crucial.  Such  a  case  is  illustrated  in 
Figure  2-18  where  transmission  losses  are  compared  for  a 
500-foot  surface  duct  at  300  Hz  (Spofford,  iy /JO) .  The 
three  predictions  were  generated  using  the  Clay/Tatro  model 
FACT  tracing  rays  in  the  duct,  and  the  NUC  surface-duct 
normal-mode  model  (Bucker,  1970)  (considered  to  be  the 
exact  solution  for  this  case) .  Sihce  the  surface  was 
assumed  to  be  flat,  and  leakage  was  not  significant  to  the 
range  shown,  the  agreement  is  perhaps  not  that  surprising. 

2.2.6  Special-Purpose  Shallow-Water  Model 

While  the  FACT  model  was  designed  for  deep-water  prob¬ 
lems,  several  of  the  FNWC  ASRAP  areas  have  bottom  depths 
less  than  1000  feet.  When  low- loss  bottom  classes  are 
specified  (i.e.  no  loss  up  to  some  apparent  critical  angle) 
a  very  large  number  of  bottom-bounce  paths  must  be  computed 
with  an  attendant  increase  in  running  time.  For  these  low 
loss  cases  an  option  has  been  provided  which  approximates 
the  contributions  of  the  surface-reflected  bottom- ref lected 
paths  with  a  rapid  computable  analytical  expression.  This 
expression  approximates  the  summation  of  all  paths  by  an 
integral  which  includes  the  surface-image  interference 
effects.  The  resulting  transmission-loss  curve  corresponds 
to  a  running  range  average  of  the  FACT  prediction. 

2.2.7  Special-Purpose  Half-Channel  Model 

Another  case  which  requires  excessive  computing  time 
for  ASRAP  forecasts  is  "half-channel”  transmission;  that  is 


2-15 


»-  V^.Tf^.1’^ ?fW5SF ‘RT'®?'*' 77***Vfl,t’'fiY?t-!r*:' 'T '^rt7W^PJ?^*?sr!5f»y r-PQftP  /-•'vt-f 


the  case  in  which  the  sound  speed  increases  montonically 
from  the  surface  to  the  bottom.  For  shallow  sources  and 
receivers  in  this  environment  (essentially  a  very  deep 
surface  duct)  considerable  computer  time  is  consumed  in 
obtaining  a  transmission  loss  curve  for  the  RSR  paths 
which  is  quite  smooth  (due  to  the  considerable  overlap  of 
tne  arrival  orders).  In  such  cases,  the  water  is  nearly 
isothermal,  and  the  key  parameters  for  RSR  propagation 
are  the  source,  receiver,  and  water  depths  and  the  frequency. 
Hence  for  the  tour  frequencies  and  three  source-receiver 
depth  combinations  used  in  ASRAP  the  FACT  RSR  intensities 
(excluding  volume  attenuation)  were  fitted  with  functions 
of  the  form  (A+BI.og  (D/1000. )  )/R  where  A  and  B  depend  upon 
the  frequency/source/receiver  depth  combination  in  question, 
D  is  the  water  depth,  and  R  the  range.  A  simple  lookup 
table  is  then  used  in  actual  runs  and  this  RSR  intensity 
is  added  to  the  direct-path  and  bottom-bounce  contributions 
which  are  computed  in  the  usual  way.  This  routine  should 
only  be  used  for  precise  ASRAP  geometries  and  frequencies 
since  it  is  not  necessarily  representative  of  other  cases. 
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Figure  2-1.  Typical  Pair  of  Paths  Experiencing  Long-Range 
Surface-Image  Interference 


Figure  2-2.  25  Hz 


Figure  2-3.  50  Hz 
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Figure  2  -S-  200  Hz 


comparison  of  FACT  ("RAY  THEORY")  with  Normal 
Modes  for  an  RSR  Convergence  Zone 
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Comparison  of  PACT  (“RAY  TRACING")  with 
Normal  Modes  for  Cusoed  Caustics 


Figure  2-14.  50  Hz 


Figure  2-15.  100  Hz 


RAY  TRACiNG 
NORMAL  MOOE 


Figure  2-13.  25  Hz 


4  /A  '\  ! 

/  (,  /|  I  /•  i;  &  \ 
\:l  tJ  -ii'.  Vi 


.80  100 


u;; 

83  ICO  '20  140 

RANGE  iw<M> 


I 


-RAY  TRACING 
NORMAL  Moot' 


A  k 
lift 

ii 1  A 


I! 


-  ,v 

*ft\  I  ■>'" 

:  I  i{ 


A 


»■  is  If  ‘ 

I .  Ii 


v*J 

80  <00  ><9 

RANGE  trie.) 


•:  1  1  •  ft- 
•  .*!■  j  \\ 

ii; 


KO  ISO  200 


Comparison  of  FACT  (“RAY  TRACING* )  with 
Normal  Modes  for  Four-Ray  Systems 
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Figure  2-16.  Ray  Paths  between  Axial  Source  and  Receiver 
for  a  Bilinear  Profile 
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Figure  2-i7,  Hay  Paths  from  a  Source  on  Axis  for  an 
Asymmetric  Quadratic  Profile 
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MATHEMATICAL  DESCRIPTION  OF  THE  FACT  MODEL 


The  previous  section  contains  the  background  physics 
and  the  rationale  for  the  approach  employed  in  FACT.  This 
section  documents  the  mathematical  implementation  of  the 
approach  and  Section  4  contains  a  detailed  flow  chart  of  the 
resulting  computer  code. 

In  this  section  some  specific  terms  are  introduced  from 
the  jargon  of  ray  tracing.  While  many  of  them  also  appear 
in  the  measurement  literature,  the  definitions  differ  some¬ 
what  and  are  repeated  here  for  conciseness.  In  ray  solutions 
for  transmission  loss  one  attempts  to  find  all  significant 
rays  or  paths  connecting  the  source  and  receiver.  The  pre¬ 
cise  rays  can,  in  principle,  be  determined  by  iterative  tar¬ 
geting  procedures,  however,  the  excessive  running  time  is 
unacceptable.  More  typically  a  number  of  rays  is  traced 
from  the  starting  point  (or  receiver)  and  the  ranges  at 
which  each  ray  crosses  the  source  depth  are  noted.  The 
source-depth  crossings  of  these  rays  are,  in  ray  jargon, 
arrivals .  If  the  ray  density  is  adequate,  the  ray  through 
the  desired  range  is  determined  by  interpolating  between 
arrivals  of  adjacent  rays  bracketing  the  range  of  interest. 
Such  interpolations  are  justified  when  rays  are  in  the  same 
family  and  of  the  same  arrival  order.  The  arrival  order 
refers  to  the  number  of  deep  cycles  the  ray  has  experienced 
(either  refracting  at  depth  or  bottom  reflecting),  within 
an  arrival  order  rays  are  subdivided  into  classes  (RR,  RSR, 
etc.).  Within  classes  additional  subdivisions  into  ray 
families  may  be  required  to  assist  the  interpolation 
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(family  3);  and  one  surface-reflected  bottom-reflected  or 
SRBR  (family  4).  In  each  family  except  the  last,  two  or 
more  rays  are  specified  within  each  layer  encompassed  by  a 
family.  The  extreme  effects  of  false  caustics  are  excluded 
by  avoiding  rays  with  horizontal  turning  points  that  are 
just  beyond  a  level  where  the  gradient  increases  in  magnitude. 
For  example,  the  ray  6^  in  Family  1  is  traced,  but  the 
slightly  steeper  ray,  0*  is  avoided.  The  final  family 

v 

(SRBR)  includes  only  the  first  ray  to  strike  both  the  surface 
and  the  bottom  and  the  steeper  of  the  following  two  rays: 

1.  The  ray  five  degrees  steeper  at  the  receiver  than 
the  first  ray; 

2.  The  ray  striking  the  bottom  at  the  effective  critical 
angle  of  the  low-frequency  bottom- ref lectivity  function. 

More  families  may  be  formed  than  those  described  in  the 

example  above  if  the  gradient,  from  one  layer  to  the  next, 

decreases  in  absolute  value  as  the  sound  speed  increases , 

thus  allowing  for  maximum-range  caustics  associated  with 

points  of  inflection  in  the  velocity  profile.  In  such  a 

case  the  rays  equivalent  to  9  and  9*  in  Figure  3-1  become 

c  c 

bounding  rays  of  the  new  families. 


When  the  rays  have  been  selected,  they  are  traced  for 
half  a  cycle  and  three  quantities  are  computed:  the  range 
to  the  first  intersection  of  the  source  depth  for  a  ray 
traced  upward  at  angle  9,  the  range  to  the  second  inter¬ 
section  for  the  same  ray,  and  the  period  of  the  ray.  From 
these  quantities  the  ranges  to  the  first  and  second  source- 
depth  intersections  for  the  downgoing  ray,  as  well  as  the 


& 
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ranges  of  all  subsequent  arrivals,  may  be  computed  using 
the  periodicity  of  the  ray  in  a  layered  environment. 


The  range  is  accumulated  as  the  ray  steps  through  each 
layer  making  use  of  the  well-known  fact  that  within  each 
layer  the  ra/  trajectory  is  the  arc  of  a  circle.  The  ex¬ 
pression  used  for  the  incremental  range,  AR,  within  a 
layer  between  depths  y,  and  y2  (with  ray  angles  Gt  and  9^ 
respectively'  i 


AR  = 


(y  -y  )*  (cos6  +cos0  ) 
2  1  2  1 

(sin0'+sin8  ) 
2  1 


(3-1) 


This  is  identical  to  the  more  traditional  expressions  and 
avoids  the  numerical  problems  associated  with  small  gradients. 
Snell's  Law,  of  course,  gives  0^  where 


cos9k  -  --  coseR 


(3-2) 


When  a  ray  turns  within  a  layer  with  gradient  g,  the  range 
increment  to  the  horizontal  turning  point  is 


AR 


c  lanh 

i  i 


(3-3) 


3.2  Smoothing  8-R  Curves 

After  the  rays  are  traced  the  arrival  ranges  of  all  the 
rays  within  a  family  are  grouped  to  form  an  arrival  order. 
The  range  to  the  arrival  as  a  function  of  the  source  angle. 
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0*  is  then  fit  with  one  of  several  functions  depending  upon 
the  family.  The  purpose  of  this  fitting  is  to  remove  false 
caustics,  obtain  the  relevant  parameters  of  true  caustics, 
and  provide  continuous  analytical  functions  for  R(9)  which 
may  subsequently  be  inverted  and  differentiated  to  obtain 
the  intensities  of  all  rays  through  all  points  of  interest. 


For  all  families  but  the  SR3R  {and  a  portion  of  families 
containing  a  cusped  caustic)  R(9)  is  fitted  with  one  of  two 
types  of  second  order  functions: 


R(9) 


(3-4) 


or  families  in  which  the  minimum  ray,  9L>  just  grazes  either 
the  surface,  the  bottom,  a  local  maximum  in  sound  speed,  or 
a  decrease  in  magnitude  of  the  gradient;  or 


R{ 9)  =  a  +  b  tan  9  +  c  tan 2 9 


(3-5) 


if  the  minimum  ray  does  not  meet  one  of  the  above  criteria. 

The  first  form  is  chosen  for  the  appropriate  families 
since  it  permits  the  presence  of  one  caustic  (R’  =  0)  and 
goes  to  the  proper  limiting  intens  'ty  (0)  as  9  ■*  0^  ( R '  -*  **») . 
This  behavior  is  characteristic  of  all  families  except 
occasionally  the  first  (shallowest).  The  second  form  may 
be  required  for  the  first  family  when  the  minimal  ray  just 
grazes  the  arrival  depth,  in  which  case  the  intensity  is 
finite  for  this  ray. 

The  coefficients  a,  b  and  c  are  determined  by  requiring 
the  function  to  pass  through  the  two  points  at  the  angular 
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extremes  of  the  family  and  through  one  of  the  following 
three  interior  points  (in  preferential  order) : 

1.  The  minimum  range  for  the  family  if  less  than  the 
range  at  both  endpoints; 

2.  The  maximum  range  if  greater  than  the  range  at  the 
endpoints; 

3.  The  second  point  in  the  family. 


For  the  last  (SRBR)  family,  R(0)  is  given  by 


a  +  b(o  eLj  ^  +c(e  -  oL) 


R  (0 )  s  {  or 


(3-6) 


d  tan  6  +  e 


The  first  expression  .is  used  for  angles  less  than  the  critical 
ray  on  the  bottom,  and  the  second  is  used  for  steeper  rays. 

The  coefficients  a  through  e  are  determined  by  requiring  R(0) 
to  pass  through  the  ranges  of  the  two  traced  rays,  continuity 
of  R' (9)  (and  hence  intensity)  at  the  critical  ray,  and  the 
proper  limiting  form  (isovelocity  medium)  for  R*  (0)  as  0  -*•  90 
degrees.  (Arrivals  of  rays  striking  the  bottom  at  more  than 
the  critical  ray  are  ignored  after  the  fourth  bottom  bounce.) 

A  subsequent  section  is  devoted  exclusively  to  the  treatment 
of  cusped  caustics  since  it  depends  upon  much  of  the  next 
section. 
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3 . 3  Computing  Ray  Intensities 


The  intensity  of  an  arrival  is  given  by; 

cos  f) 

I  =  - m-  •  (3-7) 

R  sin  o’^toL, 

where  0(R)  is  given  by  the  above  functional  forms,  O'(R) 

(the  angle  of  the  ray  at  the  receiver)  by  Snell's  Law,  and 
(3R/30)y,  by  differentiating  the  functions  for  R(0)  .  Mod¬ 
ifying  this  intensity  are  volume-absorption  and  bottom- 
reflection  losses  where  appropriate.  The  0-R  curves  for  the 
typical  arrival  order  in  the  example  of  Figure  3-1  are  shown 
in  Figure  3-2.  The  corresponding  fitted  curves  used  in  the 
program  would  differ  only  for  the  first  family  (Figure  3- 3a). 
In  this  instance  the  rigorous  ray  result  (solid  curve)  con¬ 
tains  a  true  caustic,  C,  and  a  false  caustic,  C' ,  whereas 
the  fitted  curve  (dashed)  contains  only  one  caustic.  The 
corresponding  I-R  and  total  I-R  curves  (Figures  3- 3b  and 
3- 3c)  indicate  a  significant  change  in  intensities  when  the 
false  caustic  is  removed. 

The  asymptotic  intensity  at  a  smooth  caustic  is  given  by 


1ca 


2*/*  nAi*{0)  sint//}  0  cos  0  j/, 

- - ■ — - * - “  o  ' 


l/> 


R  sin  0 ' 
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where  Ai(0)  is  the  Airy  function  evaluated  at  zero,  e  is 
the  sound  speed  at  the  source,  and  m  is  the  angular  frequency. 
The  field  in  the  vicinity  of  the  caustic  is  then 
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Ai2 (— x ) 


(3-9) 


where 


and 


I  (R) 


x 


+ 


a (R-Rc) 


(3-10) 


(3-11) 


The  sign  of  x  is  chosen  as  positive  for  points  in  the  illumin¬ 
ated  region  and  negative  for  those  in  the  shadow-zone.  The 
field  in  the  illuminated  region  is  rapidly  oscillatory  be¬ 
cause  of  the  two-ray  interference  (solid  curve  in  Figure 
3-4).  The  full  description  of  the  field  is  used  near  the 
caustic,  while  in  the  oscillatory  region  the  intensity  is 
giver,  by  the  rms  summation  of  the  two  paths  using  their 
actual  computed  intensities  (dashed  curve  in  Figure  3-4). 

The  transition  point  corresponds  to  the  first  point  (x  =  1.77) 
of  90-deqree  phase  difference  (allowing  for  the  -90  degree 
phase  shift  of  the  one  ray  which  has  touched  the  caustic). 

In  practice  the  caustic  field  is  extended  into  the  shadow 
zone  until  it  ;s  40  dB  below  the  caustic  value  (x  «  -3.5). 

The  shadow-zone  field  is  considered  to  be  propagating  at  the 
angle  of  the  caustic  ray  at  the  source  depth  (and  to  have  an 
angle  at  the  receiver  corresponding  to  the  same  ray) . 

3 . 4  Treatment  of  Cusped  Caustics 

For  cases  where  the  source  and  receiver  are  at  the 
same  depth  a  cusped  caustic  forms  on  the  zero-degree  ray  at 


each  return  to  the  receiver  depth.  Two  types  of  (smoothed) 

0-R  curves  are  possible  corresponding  to  simple  three-ray 

systems  (Figures  3-5a)  or  four-ray  systems  (Figure  3-5b) . 

In  the  latter  case  the  rays  corresponding  to  those  in  the 

one- ray  region  of  Figure  3-5a,  form  a  smooth  caustic  of  their 

own  at  (R‘ ,  6 ' ) . 

c  c 

In  both  cases  the  three-ray  region  is  treated  by  fitting 
the  branch  with  a  smooth  extremal  (connecting  points  a  and  b) 

R(C)  =  Rc  +  |2  (3-12) 

and  the  branch  joining  the  cusp  with  finite  slope  by 

R(0)  =  Rc  +  0  {  3“  1  3) 

Y 

0  and  y  ate  determined  by  the  extremals  of  the  family.  For 
the  case  of  Figure  3-5a  a  similar  fit  is  obtained  for  the 
one-ray  region  using  a  possibly  different  value  of  y. 

The  nonuniform  approximation  used  for  calculating  the 
intensity  of  the  field  near  the  cusped  caustic  is  given  by 

1  *  *  (nr)  (”k  )  [Pe (G ,  Y)f  ,  (3-14) 

where  A  is  the  wavelength,  0  is  the  ray  angle  at  the  source 
or  teceiver  (which  are  at  the  same  depth) ,  0  is  the  "cusp 
parameter",  and  Pe (0,V)  is  the  rraa  average  of  the  Pearcey 
function  Pe(X,Y)  along  its  axis  (X=0).  The  argument  Y  is 


Between  Y  =  /JtF  and  Y  »  2,  the  square  of  the  modulus  of 
Pe(0,Y)  is  used  in  tabular  form.  The  resultant  curve  for 
(Pe (6,y7) 2  is  plotted  in  Figure  3-6b. 

In  the  case  of  a  four-ray  system,  the  field  associated 
with  the  two  intersecting  branches  of  the  8-R  curve  for 
points  interior  to  the  cusp  is  computed  as  above.  For  the 
branch  forming  the  smooth  caustic  at  0^  (Figure  3-5b) ,  the 
smooth-caustic  parameters  are  computed  and  the  field  asso¬ 
ciated  with  each  portion  of  this  branch  is  estimated.  For 
the  portion  of  the  branch  corresponding  to  rays  steeper  than 
the  smooth-caustic  ray  (6  >  6^) ,  and  in  the  shadow  zone 

(R  <  R'),  the  smooth-caustic  fields  are  used.  For  the  por- 
c 

tion  of  the  branch  between  0=9’  and  8=0,  the  maximum  of 

c 

either  the  smooth-caustic  field  or  the  cusped-caustic  field 
(as  given  by  Equation  3-14)  is  used.  The  result  is  a  con¬ 
tinuous  transition  from  the  peak  field  associated  with  the 
cusp  to  the  rapidly  attenuated  field  associated  with  the 
shadow-zone  of  the  smooth  caustic. 

For  severe  cases  where  the  smooth  and  cusped  caustics 
are  separated  by  fractions  of  a  wavelength  a  more  accurate 
evaluation  of  the  field  requires  that  the  smooth-  and  cusped 
caustic  fields  be  added  on  a  phase  basis  rather  than  on  the 
rms  basis  used  here.  Such  an  approach  is  currently  beyond 
the  scope  of  FACT  and  in  such  instances  one  might  expect 
fields  slightly  stronger  (as  much  as  3  dB)  than  predicted. 

While  this  completes  the  computation  of  the  field,  the 
arrival-structure  option  requires  the  angle  and  intensity 
of  each  path  from  a  given  range.  In  the  three- ray  region. 
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it  is  assumed  that  each  of  the  three  paths  has  the  same 
intensity  and  that  the  angles  are  obtained  from  Equations 
3-12  and  3-13.  For  a  simple  cusped  caustic  the  single  ray 
carries  all  the  energy.  For  four-ray  systems  the  angles  for 
the  branch  outside  the  cusp  are  computed  as  for  any  smooth 
caustic  with  the  corresponding  partitioning  of  energy. 

3. 5  Coherent  Summation  of  Paths 

There  are  two  situations  in  which  low-frequency  coherence 
effects  are  expected  (Figure  3-7) .  In  the  first  situation 
(only  one  depth  of  interest  is  shallow) ,  the  total  inten¬ 
sity,  1^ ,  is  given  by  the  phased  sum  of  the  intensities  of 
the  two  rays  shown,  and  may  be  approximated  by 

I,  =41,  sin1  sin  9\  (3-19) 

where  I (  is  the  single-path  intensity  (assumed  equal  for 
both  paths),  and  D  the  depth  of  interest. 

When  both  depths  of  interest  are  shallow  (Figure  3-7b) , 
four  paths  must  be  phase-summed  and 

(tuD  1  sin  0  \  / WD2  sin  0\ 

— - - jsin2  y - - - J.  (3”20) 

For  large  D  and/or  u>  the  variation  in  angle,  0,  with  range 
may  produce  oscillations  which  are  inadequately  sampled  by 
the  specified  range  step.  Rewriting  the  two-path  sum  as 

I,  =  21,  ^1 


-  6  cos/2wDsin0 


“) 


(3-21) 
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we  see  that  6=1  yields  the  previous  expression  whereas  6=0 
yields  the  standard  incoherent  sum.  6  then  plays  the  role 
of  a  coherence  factor,  and  without  inferring  any  physical 
significance,  it  may  be  used  to  provide  a  smooth  transition 
from  a  fully  coherent  to  incoherent  sum  .s  the  sampling 
interval  becomes  inadequate.  For  a  given  family  the  number 
of  cycles  of  the  interference  pattern  of  frequency  f  is 
given  by 


N 

c 


2Df 

C 


(sin0 


max 


sin0  .  ) 
mm 


(3-21) 


where  0max  and  Qm^n  are  the  bounding  rays  of  the  family. 

The  number  of  range  points  (at  interval  6 R)  sampled  per  cycle 
of  interference  for  a  given  arrival  order  of  total  range-span 
AR  is  then 


N  =  AR  (3-22) 


The  smoothing,  or  coherence,  factor  6  is  then  determined  by 


0  N  i  8/3 

P 

Nd~8/3  Q/3  <  N  <  6 

“6^573  P 


6  N 


(3-23) 


The  net  effect  when  applied  to  the  2  or  4-path  cases  is  a 
smooth  transition  from  fully  coherent  to  incoherent  summa¬ 
tion  as,  for  example,  the  sampling  interval  or  the  frequency 
increases. 
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i(  '|«i  ox  i  hi. 1 1  i  on  is  doi|i. id<d  .is  I  lii>  r.iy  .iiujle  fn-eotne::  very 
shallow  of  I  ho  sour  co  (iccoivor)  hocomos  very  deep.  Tn 
such  cases  the  approximation  for  the  phase  difference  be¬ 
comes  poorer  as  does  the  assumption  of  equal  ray  amplitudes. 
Also  since  the  calculation  takes  two  or  four  families  in  an 
arrival  order  and  replaces  them  with  one  modified  family 
the  range  extent  of  the  order  is  reduced  somewhat.  To  pre¬ 
vent.  significant  errors  this  approximation  is  permitted 
only  if  the  up-  anti  down-going  arrivals  at  the  source 
(receiver)  of  the  shallowest  ray  in  the'  family  are  separ¬ 
ated  by  no  more  than  2  nautical  miles. 

i .  *.  l.ow-  Frequency  Cu  ted)  f  f.  Effects 

The  ray  equivalent  of  the  first  propaqalinq  mode 
satisfies?  t  he  phaso-intogr.il  relationship 


( 3-24) 


whei  e  t  he  •  at  eurat 


is  taken  over  a  full  ray  cycle. 


k 


Zy  sinh  (z) 


(3-25) 


and  .N s,*'  cor  responds  to  Hie  iliscrete  phase  shifts  accumulated 
by  the  ray  over  a  cycle  («  for  surface  reflections  plus  n/2 
lor  each  horizontal  reve-jal).  At  tire  relevant  sound-channel 
axis  (which  will  be  the  surface  for  half-channel  propagation) 
this  ray  has  ancle  0  .  The  intensities  of  all  rays  shallower 
than  u  are  reducoa  by  the  factor  s i n  ‘ t»x/s  i n 0Q  where  (>x  is 
the  ray  ancle  at  (he  axis. 
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3. 7  Axis-to-Axis  Propagation 


The  "period"  of  the  axial  ray  in  a  sound-speed  profile 
with  a  discontinuous  second  derivative  at  the  smooth  axis 
(i.e.  the  distance  to  F  in  Figure  2-17)  is  given  by 


'•  ■  •  (ft)’"  *  ft)'') 


(3-26) 


where  cq  is  the  axial  sound  speed,  and  c^  and  c^  are  the 
values  of  the  second  derivative  of  c(z)  just  above  and  just 
below  the  axis. 

These  derivatives  are  estimated  by  the  following  pro¬ 
cedure  (illustrated  in  Figure  3-8): 


1.  Determine  whether  or  not  the  source  and  receiver 
are  contained  in  the  two  layers  defining  the  axis  (z^) ,  and 
if  they  are,  in  the  smaller  region  defined  by  z  and  zf. 

2.  If  Zg  and  /  (z_,  z+)  exit  without  applying  axis 

corrections.  Otherwise  compute  c"  and  cM  from 

—  + 

2 (c(z*)  -  c(z  ) ) 


i.o.  fit  c ( z )  with  two  parabolas,  one  through  z  and  z  ,  the 

other  through  z  and  z  ,  both  with  their  minima  at  z  . 

♦  x  x 

3.  Determine  the  ray  which  when  traced  in  the  seg¬ 
mented  profile  has  period  P  ,  from 


tanO 

x 


(3-28) 


3-1  = 


where  g_  and  g+  are  the  gradients  above  and  below  the  axis 
respectively. 

4.  If  and  are  both  between  the  turning  points 

of  this  ray,  move  both  2  and  z  to  the  nearer  of  the  two 

b  K 

turning  points. 

In  this  procedure  the  importan*  characteristics  of  the 
axis  are  governed  by  the  two  closest  layers:  and  the  user 
should  bear  this  in  mind  when  selecting  points  for  the  input 
profile. * 

3 . 8  Sur face-bucted  Propagation 

The  transmission  loss  as  a  function  of  ranqe  R  (in  nm) 
in  the  surface  duct  (excluding  volumetric  absorption)  is 
given  by 


TL(R)  =  22  +  10Log|0  (R*A0)  +  b-R  (3-29) 


where  Ad  is  the  angle  at  the  surface  of  the  ray  just  grazing 
the  bottom  of  the  duct  and  b  contains  the  duct-leakage  and 
rough-surface  losses 


b  *  14.88x10*  f-V*  a-  7J  3,- » 


D 


(3-30) 


* / •  -»/a 

4  LJ  f  '  *  ' 

w  £k  ZD 


in  the  above  expression  t  is  the  frequency  in  kilohertz. 

K 

g  the  magnitude  of  the  be low- duct  gradient,  zR  the  duct 
d«*pth  (in  ft),  and  W  the  wave-height  factor  which  is  sea- 
state  (SS)  dependent : 


*  This  approach  to  axis-to-axis  propagation  differs  somewhat 
from  that  originally  incorporated  in  FACT  thanks  to  several 
helpful  comments  from  Mr.  C.  L.  Bartberger. 
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SS  <  3 


W  =  13.5  SS  =  3  (3-31) 

(  18  SS  >  3 

3.9  Shallow-Water  Module 

The  special-purpose  shallow-water  module  computes  the 
loss  for  cases  where  ray  computations  would  be  too  lengthy 
because  of  the  combination  of  very  shallow  water  and  a  bottom 
which  is  perfectly  reflecting  below  critical  angle.  The 
approximation  assumes  a  homogeneous  medium  and  replaces  the 
summation  of  all  paths  by  an  approximate  integral.  Surface- 
image  interference  effects  are  included  where  their  range- 
averaged  effect  leads  to  a  result  appreciably  different 
from  the  incoherent  combination  of  paths.  The  bottom  is 
treated  as  having  an  intensity  reflection  coefficient  of 
unity  up  to  the  critical  angle,  6  ,  and  a  uniform  loss  given 

v>» 

by  the  loss  at  normal  incidence,  L90,  at  steeper  angles. 

The  contributions  of  rays  steeper  than  Bc  are  included  for 
only  one  bottom- reflected  order  since  for  higher  orders 
they  are  always  negligible  compared  to  the  other  available 
paths.  The  intensity  at  any  range  then  consists  of  the  con¬ 
tributions  of  the  direct  and  bottom-reflected  paths. 

1  "  H  *  JBB  U'32! 

The  direct-path  contribution  is  computed  by  summing 

incoherently  the  direct  and  surface  reflected  paths  out  to 

the  range,  R  ,  where  the  phase  difference  between  the  two 
rms 

paths  is  -ti/2,  and  summing  coherently  beyond: 


2/R 

4  sin?/A4>\ 
R2  l2 } 


for  R  <  R 

-  rms 


for  R  >  R 


(3-33) 


where 


4tt  D  D 
=  12 

X  R 


(3-34) 


and  A <f>  =  tt/2  at  R 

rms 


8D  D 

1  2 


(3-35) 


For  bottom- bounce  paths  the  range  is  divided  into  two 
regions  R  <  R^B  and  R  >.  rqb  where 


rbb  '  2V3c 


(3-361 


is  approximately  the  period  of  the  critical  ray  in  a  homo¬ 
geneous  region  of  depth  z^.  For  R  <  RgB 


- -  •  L  *  Fl 

Ra+ ( 2Zg) 2  90 


(3-37) 


where  F1  contains  the  surface-image  inte/.ference  effects 


averaged  over  all  rays  between  9  and  n/2 

c 


p  >  B  16 


f  W/*  sir.* 

J*e  \  '  / 


sin2  telin)  ae 


(G (u j ,Uj , n/2)  —  G (u j t u j t  ) 3 


(3-38) 
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where  •  wD^sinfl/c  for  k  »  L,  2  and 
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we  obtain 


(3-44) 
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Figure  3-3.  Comparison  of  Traced  { _ )  and  Fitted  ( - ) 

Angle  and  Intensity  vs  Range  Curves  for 
Family  1  of  Figure  3-2 
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Figure  3-4.  Comparison  of  Full -Asymptotic  and  FACT  Approximation 
Treatments  of  Intensity  Near  a  Smooth  Caustic 


Figure  3- 5a.  Three- Ray 
System  of  Figure  2-8 
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Figure  3- 5b.  Four-Ray 
System  of  Figure  2-12 


«(*03tTlVC  00**1 


$-R  Curves  for  Cusped  Caustrcs 
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Figure  3- 


6.  Comparison  of  Full-Asymptotic  and  FACT- 
Approxiroation  Treatments  of  Intensity 
Near  a  Cusped  Caustic 
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Figure  3-7.  Conditions  Requiring  Coherent  Addition  at 
Low  Frequencies 
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